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Abstract We evaluate priors by the second order asymptotic behaviour of the 
corresponding estimators. Under certain regularity conditions, the risk differences 
between efficient estimators of parameters taking values in a domain D, an open 
connected subset of R'^ , are asymptotically expressed as elliptic differential forms 
depending on the asymptotic covariance matrix V. Each efficient estimator has the 
same asymptotic risk as a "local Bayes" estimate corresponding to a prior density 
p. The asymptotic decision theory of the estimators identifies the smooth prior 
densities as admissible or inadmissible, according to the existence of solutions to 
certain elliptic differential equations. The prior p is admissible if the quantity pV 
is sufficiently small near the boundary of D. We exhibit the unique admissible in- 
variant prior for V = I,D = R'^ — {0}. A detailed example is given for a normal 
mixture model. 

1 Introduction 

A parameter x takes values in a domain D, an open connected subset of i?/^. 
I use the partial differential equation symbol x rather than the usual statistical 
symbol 6 because the evaluation of asymptotic risk reduces to existence problems 
in the theory of partial differential equations. The parameter indexes a probability 
density p{yn\x) with respect to some measure say, for data y„. 

We use the KuUback-Leibler loss function 

(1) X,XGD: Ln{x,x) = j \0g\p{yn\x)/p{yn\x)]p{yn\x)dlln{yn) 

to define the risk of the estimator Xn, a function of y„ taking values in D: 

(2) Rn{Xn,x) = j Ln{Xn{yn),x)p{y„\x)dfln- 

For a prior density p, the posterior Bayes estimator x{yn,p) minimizes the posterior 
KuUback-Leibler risk 

(3) R{x\yn) = j L{x,x)p{yn\x)p{x)dx/ j p{yjx)p{x)dx. 
Define 

(4) Vn{x) = -l/j ^^log[p{yn\x)]p{yn\x)dHri{yn)- 

We will assume nVn{x) — >■ V{x), the asymptotic covariance matrix. 
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Following Bro"wn[Br79], and letting U denote the prior imiform over D, we consider 
estimators of form xijjn, U) + V^6(a;(y„, U)) for fixed decision Junctions b . D ^ D. 
Under smoothness conditions [Ha98] requiring smooth variation of the data density 
and the prior with x, the asymptotic risks for different decision functions differ only 
by terms of order n~^; therefore we define the asymptotic risk for decision function 
b, relative to the decision function 6 = corresponding to the uniform prior U, by 
the assumed limit 

(5) Rib, x) = \im^ n^[Rn{b, x) - Rn{U, x)] = , {d,{V,jbj) + \bibjVij} 
where 9, denotes the partial derivative 

For a prior with density p, the posterior bayes estimate corresponds asymptoti- 
cally to the decision function 6f = Sjlogp, and then the risk may be expressed in 
elliptic operator form 

(6) R{}f) = 2Y,.di{V,,d,^p)/^. 



It turns out that, under certain conditions of smoothness and boundedness for b 
and V, there is a risk matching prior density p for which R{b^) = R{b). Thus the 
behaviour of asymptotic risk for all smooth decisions is captured in the theory of 
elliptic differential equations, equations whose relevance to decision theory were 
first indicated in Stein[St56], but which were extensively elucidated for the normal 
location problem in Brown[Br71]. See also Strawderman and Cohen[SC71]. The 
asymptotic behavior of Bayes estimators near maximum likelihood estimators have 
been studied for loss functions of form Ln{x,x) = W{x — x) by Levit in [Le82], 
[Le83], [Le85]; in particular, he shows that the Bayes estimators form a complete 
class under certain regularity conditions. 



2 Risk matching priors 



For each decision function b we will find a risk matching prior p for which R(bP) = 
R{b). Then we need only consider decision functions of form 6^ — di \ogp 
and risks of form R{b^) = 2 di{Vijdj^Jp) / ^Jp in the asymptotic decision theory. 
This result will be proved under boundedness and smoothness assumptions using 
some standard tools from Pinsky[Pi95]. 

For the domain D with closure D , a function f is uniformly Holder continuous 
with exponent a,0<a<l, in5if 

(7) ll/llo..^= sup ^^m^<oo. 

The Holder spaces C'^'°'{D) consist of functions whose fc-th order partial derivatives 

are uniformly Holder continuous with exponent a va D. Say D' <Z<Z D ii D' 
is bounded and properly included in D. The Holder spaces C'^'"(-D) consist of 
functions that lie in C*='«(5') for each D' CC D. 

The domain D has a (7*^'" boundary dD if for each point xq G dD, there is a 
ball B centered at Xq and a 1-1 mapping tp : B ^ A C R'^, such that tp{B f) D) C 

{x e R'^,Xn > 0},tp{BndD) c{xe R'^,Xn = 0},tpG C'='°(S),^-i e C*^'"(A). 
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Condition A : D bounded, dD e C^'", Vij e C^'"(D), V positive definite in D. 
Condition A : Vij e C^'"(£)), V positive definite in D. 

We follow a standard approach which first proves results under the strong condition 
A, and then extends the results to the weak condition A by approximating D by 
an increasing sequence of bounded subdomains £)„ . 

Theorem 1: Suppose bi = dib € C^'°'{D) and assumption A holds. 
Then there exists a prior p, ^/p E C^'°'(D),p > in D, such that 

(8) R{b) = R{W) = 2 ^ . . di{Vijdj^)l^. 
Proof : 

Prom [Pi95, theorem 5.5], for some eigenvalue A, there exists 
u e C2'«(5),u > in = in dD, satisfying 

(9) 2 V . . di{Vijdju) - R{b)u = Xu. 
Since 

(10) [ Xu^ = /{2y ^^{V^J^ju) - R{b)u}u 
Jd Jd ^''1 

(11) = ^ -^J2..{ubi-2diu){ubj-2dju)Vij<0, 

it follows that A < 0. If A = 0, the corresponding eigenvector u provides the b- 
matching prior p = v? with bi = dilogp. If A < 0, from [Pi95, theorem 6.5], for 
each (f) e C^'"(L'), (j) > there exists a unique b-matching solution 
^ = ue C2'«(5), u > in £>, u = on aZ?, to the equation R{b) = R{W). 

Theorem 2: Suppose 6, = dib e C^'"(D) and assumption A holds. 
Then there exists a prior p, ,Jp e C'^'"{D),p > in D, such that 

(12) R{b) = R{V>) = 2 ^ . . diiVijdj^)/^. 
Proof: 

Select an increasing sequence of bounded domains D„ CC -D„+i,U-D„ = D. Within 
each domain, assumption A holds, so there exists a sequence of solutions 
^ = w„ e C'^'°'{Dn),Un > in £)„, such that R{b) = R{bP^),xG £>„. 

For some Xq G Di, without loss of generality set w„(a;o) = 1, all n. Note that 
any solution un,N > n is also a solution to R{b) = J?(6^"),x e -D„. A Harnack 
incquaHty[Pi95, p 124] implies that for N > n,Cn < Un{x) < Cn,x G Dn for some 
constants c„, C„. 

The Schauder interior estimate[Pi95,p86] implies, for some constant C„, for all 

N > n, 

(13) \\UN\\2,a,D„ - sup 1„ _ .„ < 

x,y€D„ ''^ 1^ y\ 

x^y 
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Thus, for each n, the sequence un is precompact in the ||wAr||2,o,D„ norm. By 
diagonahzation, there exists a subsequence, say un, that converges to u in C'^'°{D) 
for which 

(14) R{b) = R(bP) = 2^ d,{Vijdju)/u,x G D. 

Finally, we need to show that u e C^'"(-D). First note that u G C^'"(£'„), since 

(15) ||W'Ar||2,«,D„ < C„ for > n ^ ||M||2,a,C>„ < c„. 

For any D' CC D , the compact set D' is covered by UDn = D , and therefore by a 
finite subcovering, and so by a particular £)„ . Thus ||u||2,a,£)' < c„(Z)') for every 
D' CCD which implies u C C^^^iD). 

3 Explicit matching priors using the Feynman-Kac integral 

When A = in equation (9), the decision function b is the posterior bayes de- 
cision corresponding to the prior p : hi = 9,logp . We could determine the ratio 
p{x)/p{y) by the integral / hidxi over any path p connecting x and y. 

p 

When A < 0, so that h is no longer a gradient, it is plausible to attempt to find an 
approximating p by averaging these integrals over all paths between x and y. 

Consider the stochastic differential equation with initial condition Xj(0) = x : 

(16) dX,{t) = Y,.. V^^dW^it) + l{h- E,- djVij)dt 

Wc propose the Feynman Kac integral formula to specify a risk matching prior: 

(17) uix) = E[eM-Uo Ei biodX,)^piXiT))] 

where /J^ J^i ° '^■^i is the Stratonovitch stochastic integral, and T is the time 
to reach the boundary of D. I suspect that the condition A < is sufficient for 
the existence of the stochastic process and the integral when assumption A holds . 
When the formula is valid, we see that u{x) is determined as a weighted combination 

of its boundary values, with the weight at each boundary point determined by the 
path integral exp(— ^ J bidxi) over the various paths that reach that particular 

p 

boundary point. Many different priors risk matching b are available corresponding 
to different smooth assignments to the boundary values. 
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4 Decision theory for asymptotic risks 

We now apply decision theoretic classifications to the asymptotic risk formula. Our 
conclusions about the prior p will depend only on the domain D and the asymp- 
totic variance V . From now on we will drop the term asymptotic. We consider a 
particular variance V and the set of risks, real valued functions on the domain D, 
corresponding to priors satisfying 

(18) Assumption B : p G C'^'"{D), Vij € C^'"(L»);j9, V > in D. 
The risks will be written R{p) = RibP). 

The posterior bayes decision is locally Bayes: for any alternative decision 

+ V where v e C^'"(£>) and u = in D — D' , some D' CC D, integration by 
parts shows 

(19) / [R{lf + v)-R{hP)]p>Q. 

Jd 

Theorem 3: Under assumption B, with D bounded, the following conditions are 
equivalent: 

iP is Bayes: there exists no p* ^ p with J {R{p*) — R{p))p < 0. 
If is Admissible: there exists no p* ^ p with R{p*) < R{p),x £ D. 

is Unique Risk: there exists no p* p with R{p*) = R{p), x G D. 
W is Brown: no non-trivial positive h solves di{pVijdjh) = 0. 
Proof: 

Bayes implies Admissible because u* ^ u violating Admissible also violates Bayes. 
Admissible implies Unique Risk because u* ^ u violating Unique Risk also vi- 
olates Admissible. Brown[Br71,1.3.9] and Unique Risk are equivalent, because 
R{py/h) = Rip) if and only if J2^J d,ipV^,d,h) = 0. 

It only remains to show that failure of Bayes implies failure of Unique Risk. 
Without loss of generality, assume p = U is uniform, so that R{p) = 0. If p is not 
Bayes, there exists p* ^ p with 

(20) > j Rip*) = I E., + E., d^iVijb*) 
Since p* ^ p and the middle integral is positive, then 

(21) / ^..5i(y,,-6*)=C<0. 

JD — ^3 

Let \dD\ be the lebesguc measure of the boundary dD when D is smooth and 
bounded, let r denote the outward pointing normals on the boundary, and note 
that 

(22) / ^{r,6*y.,} = I ^^(^^^^*) = ^ < 0- 
JdD Jd '■I 

Applying Theorem 6.31 from [GT97], for D„ CC D, since C < 0, 'Etjinr^V^j} > 0, 
there exists a solution Pn = & C^'"(-D„) to the oblique derivative problem 

(23) Ripn) = in D„,w„ = 2 * C"Y{TidiuVij}/\dDn\ in dD^ 

ij 

SO that Jjj diiVijdiPn/Pn) = C < 0. 
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Repeating the compactness argument of theorem 2 on £)„ CC £)„+i,UD„ = D, 
there exists po e C'^'°'{D) with 

(24) Ripo) = 0, / V . . a,(F,,9iPo/Po) =C<0. 

The first condition states that po and p = U have the same risk, and the second 
condition guarantees that po ^ p = U, so that the unique risk condition fails, as 
required. 



5 When is pV Bayes on R'' ? 

Brown's condition shows that the locally Bayes estimate is Baycs or not de- 
pending only on the product pV. For example, is Baycs with V, if and only if 
is Baycs with pV. We will therefore rephrase the admissibility question in terms of 
the product pV: the prior-scaled covariance matrix pV is Bayes on D if and only 
if there is no non-trivial solution to Brown's equation. 

Theorem 4. Let D = , pV e C^'°'{D),r = \x\,x = rs where s ranges over the 
surface S of the unit sphere. Define W{R,s) = [J^ ^V~^r^~'^dr]~^. Suppose that, 
uniformly over s Cz S, 

(25) lim W{R,s) ^ W{s), lim W{s)W~\R, s)W{s) = W{s). 

Then pV is Bayes on only if 



(26) 



/ y^^SiSjWij{s)ds = 0. 
Jses 



Proof: 

Let Q = {Qi, 1 < i < d} G C^'"{D) be a test function with relative risk 
(27) R{bP + Q)- R{b) = J2^di{pVijQj) + ^Q.QjpVij}. 

ij 

Prom theorem 2 for every Q there exists a prior q with R{IP + Q) = R{b'^). Thus 
pV is Bayes if and only if 



(28) 



{R{lf + Q)- R{bP)}p = [ V{5i(pFi,Qj) + iQiQjpVij} > 



for every test Q where the integral is defined. Equivalently, with the test (pV) ^Q, 

(29) 

t{Q) = I {R{W + {pV)-'Q) - R{bP)}p = f {J^diQi + I Eij QiQjV^'/p} > 
Jd Jd ^ 

for every test Q where the integral is defined. 

The possible negative term ^ • diQi is determined by values in the neighbour- 
hood of infinity, so pV being Baycs is determined by the behaviour of pV near the 
infinite boundary. In particular if two functions pV are identical outside a compact 
subset of D, they have the same admissibility classification. 
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We therefore consider a test function that is zero inside the unit sphere: 
(30) Q{rs) = g{ry-'^q{s), q{s) = -Ws 

where g is twice differentiable, g{r) = for < r < 1,0 < g{r) < 1 for 1 < r < 
2, g{r) = 1 for r > 2. 

(31) 

Let t{Q, R)= [ {V diQi + \ Y.ij QiQjV^T^IvV-^drds = J^^^ I{R, s)ds. 

J\x\<R i 

Consider the contribution I{s, R) to the test integral for a particular s: 

(32) I{s,R) = VA-,Q,(i?s)i?''-i+ / {lE^,Q^QJVij'/p}r'^~'dr 

i Jo 

(33) = ^s,g,(s) + ^5^ftg, rff(r)V-%Vpdr 

i ij 

(34) < -^s,s.T^,. + ^^s,s,-W^,feW„W,-^i(i?,5) 

ij ijkl 

(35) — >■ — - SiSjWij uniformly in s as J? — > cxd 

ij 

Thus 

(36) t{Q) = Urn / /(s, R)ds < 
unless 

(37) / y2siSjWij{s)ds = 0. 

Jses 

which shows that the condition in the theorem is necessary for p to be Bayes. 

Failure of the condition in the theorem allows construction of an explicit test 
function for showing p to be not Bayes. 1 suspect that the weaker condition 
]imii^^{J2ij ^i^AIi^ ^V~^r^~'^dr]~^ds = is also necessary. It may be that the 
condition is also sufficient. A similar condition for the recurrence of diffusion pro- 
cesses is given in [Ic78]. 

Brown[Br71] studies the admissibility of estimates for the normal location problem 
in d dimensions in which it is assumed that the data x are gaussian with unknown 
mean and identity covariance matrix. He shows that an estimate corresponding to 
the marginal density of the data p{x) is admissible if Pi = dilogp = -£^\ogp is 
bounded and if 

(38) / [/ pds]-^r^-'^dr = oo. 

Ji Jses 
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Brown, Theorem 6.4.4, also shows that an estimate corresponding to the marginal 
density p{x) is admissible only if 

(39) / [/ p-^r^-'^dr]-'^ds = 0. 

J ses Ji 

The asymptotic version requires data X„ ~ N{9,I/n), with n — oo. Theorem 4 
implies (40): a prior p is Bayes only if 

(40) / [/ p-^r^-'^dr]-^ds = 

Jses Ji 

or equivalently, almost everywhere on S 

/oo 
p-^r^-ddr = oo. 

If the prior density p is expressed as a density p on the polar co-ordinates x — rs, 
the condition simplifies to P~^{r, s)dr = oo almost everywhere on S. See Straw- 
derman and Cohen[SC71], theorem 4.4.1. For example, the prior corresponding to 
r" being uniformly distributed is Bayes in every dimension for a < 2 but not Bayes 
for a > 2. 



7 When is V Bayes on bounded D? 

Let D be a bounded domain with boundary in C^'" . Let iy{s),s E dD denote 
the outward pointing normal at a point s on the boundary of D, assumed defined 
almost everywhere in ds, lebesgue measure on the boundary. It will be assumed 
that, for almost all s E dD, the inward pointing normal{s — uiy{s), < u < s} lies 
in D for e small enough. 

Theorem 5. The covariance matrix pV is Bayes only if 

(42) lim / ViVil [ ^V,~Hs-uiy)du]^'^ds = 0. 

This is proved similarly to theorem 4, using test functions that are constant on the 
inward pointing normal segments. 

It may happen that, for each s, the normal vector ^{s) ad, s is the limit of some 
eigenvector of pV as x G D -^sG dD, (the normal eigenvector case) in which case 
the condition in the theorem simplifies to 

(43) / ^^^VjJ^{s — uv)du = oo almost all s. 

Jo P 

We will say that the integral condition fails at s if /J ^^^V~^'^ {s — uv)du < oo. The 
theorem now states that pV is admissible only if the integral condition fails on a 
set of measure 0. 
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The admissible pV are those where pV — > fast enough near the boundary. If pV 
is inadmissible, we can render it admissible by attenuating p near the boundary. In 
the normal eigenvector case, let consist of those points x £ D within e of the 
boundary, and suppose that each such point is closest to a unique boundary point 
s{x). Each such point may be written x = s(x) — u{x)i'{x) for some u. 
Let a be an attenuation factor defined at each point in D by: 

a{x) = 1 ioY X <E D — D^, 

a{x) = 1 for X G D^, integral condition holds at x(s), 

a{x) = mm[l, 1 — (1 — ^ ^ j^^^ integral condition fails at x{s). 



The proposed attenuation factor will be 1 except near boundary points s where the 
integral condition fails, where it will approach zero. With the prior ap, the integral 
condition becomes 




9 One dimension 

For the one dimensional parameter x on D = (a, b) with variance V, Brown's 
condition implies that pV is admissible if and only if 

(45) ^ iVp)-' = j (Vp)-^ = 00. 

Since a smooth monotone transformation renders pV equal to 1 on _D = (a, &), an 
equivalent result is that there exists a non-zero differentiable test function w on D 
such that f^{w' + ^ui^) < if and only if either a or 6 are finite. 

Jeffreys' density J = is admissible on D = (a, b) if and only if 

(46) j j = 

which means that Jeffreys must be "improper" in both tails to be admissible. I take 
a certain delight in this impropriety, because although "improper" priors abound 
in decision theory and in Bayesian analysis, they remain objects of suspicion. See 
for example the excellent review in [KW96]. However, in decision theory, the prior 
appears only when multiplied with a loss function, which may be arbitrarily scaled, 
so improper priors form a natural part of the range of procedures we need to study. 
Asymptotically, the prior appears only as a product with the covariance matrix 
in admissibility questions, and again it makes no sense to constrain priors to be 
improper. In the Jeffreys' case, the admissibility of the product pV requires that 
Jeffreys be improper in the tails. 
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The pearson correlation coefficient computed for n bivariate normal observations 
with true correlation p has asymptotic variance l/(f — /C^)^ . 
Thus a prior p on D = (—1, 1) is admissible if and only if 

(47) j ^{{l- p-'fp)-^ = f ((1 - p2)2p)-i = oo. 

For priors of form p=l/(l — p^)", the prior p is admissible if and only if a < 1 . 
Thus if we wish to skirt the edge of inadmissibility, we might use p = l/{l — p^) . 

10 Invariant admissible priors 



Since the KuUback-Lciblcr loss fuiic;tion docs not change under smooth transfor- 
mation of the parameter space, differences between the asymptotic risk functions 
for two priors are also invariant under such transformations. We are free to trans- 
form to a convenient p, V. D in deciding admissibility problcnns. If a transforma- 
tion T takes p, V, D into say T{p), T{V), T{D), then the admissibility of PV in D 
equals the admissibility of T{p)V{p) in T{D). A prior p is relatively invariant if 
p{Tx)J — cp{x), X E D, where J is the Jacobian of the the transformation x Tx 
and when T is one to one D ^ D such that TV{Tx) = CV{x). 

For example ii D = — 0,V = I, arbitrary rotations and scalings leave D in- 
variant, and change the covariance by a constant, so the only invariant priors are 
of form J) = r",r^ = xf. From Brown's condition, the prior p is admissible if 

(48) / r^-'^-"dr = 
Jo 

which occurs only when a = d — 2. In this case there is a single admissible invariant 
prior p = r**"^. This prior, discussed in [Br71] and [SC71], corresponds to being 
uniform over D. 

If Z) = {a;|i?i < r = |a;| < i?2}, the invariant transformations are rotations, which 
require that an invariant p depends only on r. Admissibility requires 

f 1 f^^ 1 

(49) / -dr= -dr = oo. 

JriP J P 

Admissibility is achieved by p{x) = miny^dD \x — y\- 

Although invariancc considerations no longer always apply, the above solution can 
be extended to general bounded D with V = I, namely p{x) = mixiy^gD \x — y\. 
For general D,V, define |a; — y| as the path length between points x,y G D in 
the metric d{x,y) = {x — yyv~^{x — y). Then D is bounded in this metric if all 
paths have finite length, and we again define p{x) = min^goij |a; — y|. We offer this 
merely as a suggestion for an admissible prior that flirts with inadmissiblity near 
the boundaries, and is consistent under transformations of the data. 
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11 A mixture model 

Suppose that t/„ is a sample of size n from the normal mixture 

(50) Y = Z+{l-B{q))x,+B{q)x2 

where Z ~ iV(0, l),B{q) ~ Bernoulli with mean q, xi > 0,X2 > 0,q — xi+x2 ' 
The parameter x = (a;i, X2) lies in the domain D = {x^ > 0,X2 > 0}. 
The density of a single observation y is 

(51) f{y) = {X2<l){y + xi) + xi4>{y - X2)}/(xi + X2). 

The asymptotic variance V is the inverse of the information matrix of expected 
values of products of the score functions: 

.52) I 1 ^ 4){y-X2) {y + xx)x2(i){y + xi) 

^ X1+X2 {Xl+X2)f {x\-\-X2)f 

1 , (Phl + xi) {y ~ X2)xi(j){y - X2) 



^^^^ X1+X2 ' {Xi+X2)f ' {Xl + X2)f 

(54) Lij = j kljfdy 

(55) y = L-i 

Asymptotic admissibility for the prior p is determined by behaviour of Vp near the 
boundaries. Let xi = rsi, X2 = rs2, sf + S2 = 1- 

(56) a;i ^- : ill (exp(a;i) - 1 - x^/xj, L^2 0, L22 0, 

(57) a;2^0:-£/22 {e-xp{x\) - I - xl)/x\,Li2 ^ Q,Lii ^ 

(58) 00: ill -)■ S2/(Sl +S2),-f'12 ^ 0,L22 si/(si +52)- 

At the boundary xi = the normal is an eigenvector at all points (0, X2), and the 
integral condition for admissibility for that boundary is ^Ludxi = 00 almost 
all X2 which reduces to Jq ^dxi = 00 almost all X2- Similarly, the condition for 

admissibility on the boundary 0:2 = is ^dx2 = 00 almost all Xi. 

For the infinite "boundary" r ^ 00, the integral condition for admissibility is 

(59) lim / SiSj[riv^\rs)r-'dr]-'ds = 0, 

where S is the intersection of the boundary of the unit circle and the upper right 
quadrant. Using the behavior of L as r ^ 00, this condition becomes 

(60) / — z — rdr — >■ 00 almost all s G S. 

Ji rp[sr) 

Choosing a prior p to make pV admissible requires that 

(61) / -dxi = / -dx2, / — dr = 00. 
Jo P Jo P Ji 

Roughly, we need that p be of order xi near .ti = , of order X2 near a:;2 = ,and 

of order log{r) near r = 00 . For example, p = (^xi+x2)^ ^^^^ ^ ^^^^ 

many other priors with the correct behavior near the boundary. The uniform is 

inadmissible because it fails at a;i = and a;2 = 0. 



Confidence ellipsoids for x1,x2 



X 




x1 



Figure 1. Confidence Ellipses 

The plot of confidence ellipses when 1000 points are sampled from the mixture 
model shows how the boundaries affect asymptotic admissibility. For the bound- 
aries xl = and x2 = 0, the asymptotic variances orthogonal to the boundary 
in fact approach a positive limit; thus the integral of the inverse variances up to 
the boundary is positive rather than infinite, and the uniform density is therefore 
inadmissible. For the boundary at infinity, the variances are bounded away from 
infinity, so the integral of the inverse variance is infinite, and this boundary is 
admissible for a uniform prior. 
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12 A prior beating the uniform 

It is of interest to exhibit a prior with asymptotic risk everywhere smaller than 
an inadmissible prior such as the the uniform in this problem. Brown's condition 
exhibits a prior satisfying ^[V^) = 0. The asymptotic risk of p, relative to the 

uniform, is ^5(gf)'^§f/p • There arc many solutions to the elliptic differential 
equation, depending on boundary values of p. The solutions are not necessarily 
admissible. 

We have computed solutions for the discrete approximation where Xi , X2 each lie 
in the grid 0.1, 0.2,. . . 10. The boundary values for p are p = min( (^^^^^^^2 , XiX2)- 
We set these values so that p will satisfy the conditions for admissibility at the dif- 
ferent boundaries. The following prior is obtained by using a relaxation method to 
solve the finite difference form of the differential equation; at the solution, the finite 
difference expressions are everywhere less than .01. A similar prior was developed 
in [Em02]. 

It will be noted that the prior density approaches zero at the lower and left 
boundary, but not at the other two boundaries, as required by the admissibility 
conditions. 

The risk gains against the uniform are everywhere positive ( as required by the 
theory), but are far greater near the low xl and x2 boundaries. This is to be ex- 
pected, because the prior is made admissible by changes near the boundaries, so 
that larger improvements in the risk should occur there. 
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Figure 3. Risk Gains compared to Uniform 
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